arXiv: 1503.04468vl [math.OC] 15 Mar 2015 


Ordered Median Hub Location Problems with Capacity 

Constraints 

J. Puerto^ A.B. Ramos“ A.M. Rodnguez-Chia^ M.C. Sanchez-GiP 

“ Facultad de Matemiticas. Universidad de Sevilla, Spain. 

^ Facultad de Ciencias. Universidad de Cadiz, Spain. 

March 17, 2015 


Abstract 

The Single Allocation Ordered Median Hub Location problem is a recent hub model 
introduced in |36| that provides a unifying analysis of a wide class of hub location mod¬ 
els. In this paper, we deal with the capacitated version of this problem, presenting two 
formulations as well as some preprocessing phases for fixing variables. In addition, a 
strengthening of one of these formulations is also studied through the use of some fami¬ 
lies of valid inequalities. A battery of test problems with data taken from the AP library 
are solved where it is shown that the running times have been significantly reduced with 
the improvements presented in the paper. 


1 Introduction 

Network design problems are among the most interesting models in combinatorial optimiza¬ 
tion. In the last years researchers have devoted a lot of attention to a particular member 
within this family, namely the hub location problem, that combines network design and lo¬ 
cation aspects of supply chain models, see the surveys mmm- The main advantage of using 
hubs in distribution problems is that they allow to consolidate shipments in order to reduce 
transportation costs by applying economies of scale; which are naturally incorporated to the 
models through discount factors. Hub location problems have been studied from different per- 
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spectives giving rise to a number of papers considering different criteria to be optimized: the 
minimization of the overall transportation cost (sum) (see [ini 01201 [271 |29l|30l EH), the min¬ 
imization of the largest transportation cost or the coverage cost ( 0 0 El [251120 El IIQIIII] ), 
et cetera. 

Apart from the choice of the optimization criterion, another crucial aspect in the literature 
on hub location, and in general on any location problem, is the assumption of capacity con¬ 
straints. One can recognize that although this assumption implies more realistic models, the 
difficulty to solve them also increases in orders of magnitude with respect to their uncapaci¬ 
tated counterpart, fn many cases new formulations are needed and a more specialized analysis 
is often required to solve even smaller sizes than those previously addressed for the uncapaci¬ 
tated versions of the problems. For this reason, capacitated versions of hub location problems 
have attracted the interest of locators in the last years, see [210[II1[I0[I11[I0[I0E0[29]. 
fn the same line, we also mention some other references related with congestion at hubs, as 
congestion acts as a limit on capacity, see [13 [10 Eg. 

An interesting version of hub location model is the Capacitated Hub Location Problem 
with Single Allocation (CSA-HLP), see [111 1131 EH] , fn this context, single allocation means 
that incoming and outgoing flow of each site must be shipped via the same hub. fn contrast to 
single allocation models, where binary variables are required in the allocation phase, multiple 
allocation allows different delivery patterns which in turns implies the use of continuous 
variables simplifying the problems. The CSA-HLP model incorporates capacity constraints 
on the incoming flow at the hubs coming from origin sites or even simpler, on the number of 
non-hub nodes assigned to each hub. The inclusion of capacity constraints make these models 
challenging from a theoretical point of view. Regarding its applicability we cite one example 
described in Ernst at al. |20j based on a postal delivery application, where a set of n postal 
districts (corresponding to postcode districts represented by nodes) exchange daily mail. The 
mail between all the pairs of nodes must be routed via one or at most two mail consolidation 
centers (hubs), fn order to meet time constraints, only a limited amount of mail could be 
sorted at each sorting center (mail is just sorted once, when it arrives to the first hub from 
origin sites). Hence, there are capacity restrictions on the incoming mail that must be sorted. 
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The problem requires to choose the number and location of hubs, as well as to determine the 
distribution pattern of the mail. 

The CSA-HLP has received less attention in the literature than its uncapacitated counter¬ 
part. Campbell presented the first integer Mathematical Programming formulation for the 
Capacitated Hub Location Problem. This formulation was strenghthened by Skorin-Kapov 
et al. |39) . Ernst and Krishnamoorthy |20) . proposed a new model involving three-index 
continuous variables and developed a solution approach based on Simulated Annealing where 
the bounds obtained are embedded in a branch-and-bound procedure devised for solving the 
problem optimally. Recently, Correia et al. m have shown that this formulation may be 
incomplete and an additional set of inequalities is proposed to assure the validity of the model 
in all situations. A new formulation using only two indices variables was proposed by Labbe 
et al. m, where a polyhedral analysis and new valid inequalities were addressed. Although 
this formulation has only a quadratic number of variables, it has an exponential number of 
constraints, and to solve it the authors developed a branch-and-cut algorithm based on their 
polyhedral analysis. Contreras et al. m presented for the same problem a Lagrangean re¬ 
laxation enhanced with reduction tests that allows the computation of tight upper and lower 
bounds for a large set of instances. 

In two recent papers, IMIET], a new model of hub location, namely the Single Allocation 
Ordered Median hub location problem (SA-OMHLP), has been introduced and analyzed. 
This problem can be seen as a powerful tool from a modeling point of view since it allows a 
common framework to represent many of the previously considered criteria in the literature 
of hub location. Moreover, this approach is a natural way to represent the differentiation of 
the roles played by the different parties (origins, hubs and destinations) in logistics networks 
[m EH ESI |32]. This model does not assume, in advance, any particular structure on the 
network mm)- Instead of that this structure is derived from the choice of the parameters 
dehning the objective function. Apart from the above mentioned characteristics, ordered 
median objectives are also useful to obtain robust solutions in hub problems by applying k- 
centrum, trimmed-mean or anti-trimmed-mean criteria. It is worth mentioning that although 
it is called single allocation, its meaning slightly differs from the classical interpretation in 
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hub location where each site is allocated to just one hub and all the incoming and outgoing 
flow to-from this site is shipped via the same link (the one joining this site and its allocated 
hub). In this model, single allocation means that all the outgoing flow is delivered through 
the same hub, but the incoming flow can come from different hubs. Actually, this is a mixed 
model and basically the same situation described above, about postal deliveries, naturally fits 
in this framework assuming that letters from the same origin should be sorted, with respect 
to their destinations, in the same place and from there they are delivered via their cheapest 
routes. Observe that in this scheme it is also natural that incoming flow in a final destination 
comes from different hubs. 

The SA-OMHLP distinguishes among segmented origin-destination deliveries giving dif¬ 
ferent scaling factors to the origin-hub, hub-hub and hub-destination links. The cost of each 
origin-hrst hub link is scaled by a factor that depends on the position of this cost in the 
ordered sequence of costs from each origin to its corresponding first hub OSSIES]. More¬ 
over, the overall interhub cost and hub-destination cost are multiplied by other economy of 
scale factors. The goal is to minimize the overall shipping cost under the above weighting 
scheme. The reader may note that the first type of scaling factors mentioned above adds a 
“sorting” problem to the underlying hub location model, making its formulation and solution 
much more challenging. This model and two different formulations were introduced in |35| 
while a specialized B&B&Cut algorithm was developed in |36[I38| . None of those formulations 
could handle capacities since the computation burden of the problems were highly demanding. 
Thus, the SA-OMHLP with capacity constraints, i.e. Capacitated Single Allocation Ordered 
Median hub location problem (CSA-OMHLP) is currently an open line of research for further 
analysis. 

In this paper, we analyze in depth the CSA-OMHLP trying to obtain a better knowl¬ 
edge and alternative ways to solve it. Thus, the contributions of this paper are threefold. 
First, it combines for the hrst time three challenging elements in location analysis: hub facil¬ 
ities, capacities and ordered median objectives; proposing a promising IP formulation which 
remarkably reduces the number of decision variables. Second, this paper strengthens that 
formulation with variable fixing and some families of valid inequalities that have not been 
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considered before. Finally, despite the difficulty of considering simultaneously capacitated 
models, hubs and ordering, the techniques proposed in this paper allow to solve instances of 
similar sizes to those already considered in the literature for simpler models (uncapacitated 
and multiple allocation [22]). 

The paper is organized as follows: in Section[2]we will provide, first, a MIP formulation for 
the capacitated version of the problem extending the one in |36| and then another formulation 
in the spirit of |3Z| where the number of variables has been considerably reduced with respect 
to the previous one. Section [3] strengthens the latter formulations with variable fixing and 
several new families of valid inequalities. In Section Ul the effectiveness of the proposed 
methodology is shown with an extensive computational experience comparing the performance 
of the two formulations and the strengthening proposed along the paper. Finally, the paper 
ends with some conclusions. 

2 Model and MIP formulations 

The goal of this paper is to analyze the CSA-OMHLP. For this reason, we elaborate from 
the most promising formulations of the non-capacitated version of that problem, namely the 
so called radius (covering) formulations, see |361137] . In order to be self-contained and for 
the sake of readability, we include next a concise description of these formulations in their 
application to the capacitated problem. 

Let A = {1,... ,n} be a set of n client sites, where each site is collecting or gathering 
some commodity that must be sent to the remaining ones. It is assumed, without loss of 
generality, that the set of candidate sites for establishing hubs is also A. Let Wjm > 0 be the 
amount of commodity to be supplied from the j-th to the m-th site for all j, m ^ A, and 
let Wj = Ylm£A "^jm be the total amount of commodity to be sent from the j-th site. Let 
Cjm P 0 denote the unit cost of sending commodity from site j to site m (not necessarily 
satisfying the triangular inequality). It is assumed free self-service, i.e., Cjj = 0, Vj G A. 

Let p < n be the number of hubs to be located and let hj be the capacity of a hub located 
at site j, with j ^ A. A solution for the problem is a set of sites X C A with \X\ = p and 
enough capacity to cover the flow coming from the sites; plus a set of links connecting pairs 
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(flow patterns) of sites j, m for all j, m (z A. Moreover, it is assumed that the flow pattern 
between each pair of sites traverses at least one and no more than two hubs from X. 

As it was mentioned in Section [H the main advantage of using hubs is to reduce costs 
by applying economies of scale to consolidated flows in some part of the network. In this 
model the transportation cost is decoupled into the three differentiated possible links: origin 
site-first hub, hub-to-hub, and hubs-final destination. These transportation costs are scaled 
in a different way. The model weights origin site-first hub transportation costs by using 
parameters A = (Ai,..., Xn), with A* > 0 Vi G A, depending on their ordered rank values. 
This is, let Cjk be the cost of the overall flow sent from the origin site j if it were delivered 
via the first hub k, i.e. Cjk ■= Cj^Wj, j,k G A. Next, if Cjk were ranked in the i-th position 
among all these costs, then this term would be scaled by A* in the objective function. For 
the two remaining links there is a compensation factor 0 < /r < 1 for the deliveries between 
hubs, and another one 0<(5<1,^<5, for the deliveries between hubs and final destination 
sites. These parameters may imply that, even in the case where the costs satisfy the triangle 
inequality, using a second hub results in a cheaper connection than going directly from the 
first hub to the final destination. Actually, it represents the application of the economy of 
scale by the consolidation of flow in the hubs. 

In the following we present a first valid formulation of the CSA-OMHLP, based on covering 
variables (the reader is referred to |36l I37| for further details). Sorting the different delivery 
costs values (cjk) for j, A; G A, in increasing order, we get the ordered cost sequence: 


C(l) 


0 < C( 2 ) < • • • < C(G) := max {cjk}- 

l<j,k<n 


where G is the number of different elements of the above cost sequence. For convenience we 
consider C(o) •= 0- 

For i G A and h = 1,... ,G, we define the following set of covering variables. 


^ih ■ ^ 


if the i-th smallest allocation cost is at least C(/j), 
otherwise. 


( 1 ) 


Clearly, the i-th smallest allocation cost is equal to C(/j) if and only if Uih = 1 and = 0. 
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In addition, this formulation uses two more sets of variables: 


Xjk 


Skim 


1, if the commodity sent from origin site j goes first to the hub k, 

( 2 ) 

0, otherwise. 

flow that goes through a first hub k and a second hub I with destination m, 


with j,k,i,m G A. Since we assume that any origin is allocated to itself if it is a hub, the 
above definition implies that site k is opened as a hub if the corresponding variable Xkk takes 
the value 1. 


The formulation of the model is: 

G 

{Fu) min EE Ai(c(ft,) - C(h_i))uih + EEE {pCkl + dcem)skem (3) 

ieAh=2 k£Al&Am£A 

s.t. = Vj G A (4) 

k£A 

y^^Xjk < nxkk, '^k £ A (5) 

jeA 

'^Xkk=P ( 6 ) 

keA 

^klm = E WjraXjki ^k,m€A (7) 

i&A jeA 

Skim E ^ ^ S^mm) V/c, £,171 £ A, £ ^ 171 (8) 

jeA 

EE ■^kim — ^kk E Wj, yk£ A ( 9 ) 

leAm^A jeA 

EE ^kim — E Wj, y££A (10) 

k^AraGA j^A 

jeA 

n 

= ^ ^ Xjk, yh = l,...,G (12) 

ieA jeA k=i 

Cjk>C(h) 

Uih>Ui-i,h, yi£A\{l},h = l,...,G (13) 

Uih,Xjk£ {0,1}, Skim >0, yi,j,k,£,m £ A, h = l,...,G (14) 


The objective function ([3]) accounts for the weighted sum of the three components of the 
shipping cost, namely origin-first hub, hub-hub and hub-destination. The origin-hub costs are 
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G 

accounted after assigning the lambda parameters, i.e. ^ ^ In addition, 

zGA h=2 

the second and third blocks of delivery costs, i.e. the hub-hub and hub-destination cost, 
scaled with the /r and 5 parameters respectively, can be stated as: ’Yhk&A Yli&A YlmeAiP^ke + 

The constraints of the model are described in |36) with the only exception of (llip . that 
is the capacity constraint on the hubs, in spite of that and for the sake of completeness we 
include below a brief description of them. Constraints (|4| ensure that the flow from the origin 
site j is associated with a unique first hub. Constraints ([5|) ensure that any origin only can be 
allocated to an open hub. Constraint ([^ fixes the number of hubs to be located. Constraints 
© are flow conservation constraints, such that the flow that enters any hub k with final 
destination m is the same that the flow that leaves hub k with destination m. Constraints 
([8|) ensure that if the final destination site is a hub, then the flow goes at most through 
one additional hub. These constraints are redundant whenever the cost structure satisfies the 
triangular inequality, however they are useful in reducing solution times (see m)- Constraints 
([9|) and (IlOp establish again that the intermediate nodes in any origin-destination path should 
be open hubs. Constraints (lllh establish the capacity constraints of the hubs. Observe that 
this family of constraints make redundant the family ([5]), but we have kept it because it 
reduces the computational times. Constraints (1121) link sorting and covering variables. They 
state that the number of allocations with a cost at least C(/j) must be equal to the number 
of sites that support shipping costs to the first hub greater than or equal to Finally, 

constraints m are a group of sorting conditions on the variables Uih- 

The reader may note that this formulation is a natural extension for the capacitated 
version of the radius formulation already considered for the uncapacitated ordered median 
hub location problem in IMEZI. However, although this formulation is enough to specify 
the CSA-OMHLP, we have found that for solving medium sized problems it produces very 
large MIP models, which are difficult to solve with standard MIP solvers (CPLEX, XPRESS; 
Gurobi...). Therefore, some alternatives should be investigated. 

One way to improve the performance of the above formulation is to take advantage of 
some features of that model to reduce the number of variables. In this case, one can sue- 
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ceed reducing the number of u variables. The logic of the above formulation can be further 
strengthen for important particular cases of the discrete ordered median hub location prob¬ 
lem. In the following, we show a reformulation that is based on taking advantage of sequences 
of repetitions in the A-vector. (See [331 l37l l38] for similar reformulations applied to other 
location problems.) 

One can realize that for A-vectors with sequences of repetitions -i.e. the center, /c-centrum, 
trimmed means or median among others-, many variables used in formulation Fu are not 
necessary (since they are multiplied by zero in the objective function), and some others can 
be glued together (since they have the same coefficient in the objective function). Moreover, 
under the assumption of the free self-service, and that any origin is allocated to itself if it 
is a hub, we conclude that the p smallest transportation costs from the origin to the first 
hubs are 0, i.e. the first p components of the A-vector are multiplied by 0. Therefore, in 
order to simplify the problem one can disregard the p first components of the A-vector. Let 

A — (Ai,..., Afj—p) .— (Ap^i,..., Afj). 

in order to give a formulation for the CSA-OMHLP taking advantage of these facts, we 
need to introduce some additional notation. Let I be the number of blocks of consecutive 
equal non-null elements in A and define the vectors: 

1 . 7 = ( 71 ,... , 7 /), being ji, i = 1 ,..., / the value of the elements in the Lth block of 
repeated elements in A. 

2. a = (ai,..., ai, 0 : 7 + 1 ), being Oj with i = 1 ,..., /, the number of zero entries between 
the {i — l)-th and f-th blocks of positive elements in A and 07+1 the number of zeros, 
if any, after the J-th block of non-null elements in A. For notation purposes we define 
cxq = 0 . 

3. /3 = {Pi,..., /3i), being Pi, i = 1,...,/ the number of elements in the z-th block of 
non-null elements in A. For the sake of compactness, let Pq = Pi+i = 0. 

Next, let denote oii = Pi — Z]}=i/^i recall that Wi = Moreover, for 
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all i = 1 ,...,/ and h = 1,... ,G, let us define the following set of decision variables: 


Uih 


Vih 


1 , if the {p + ai + I3^_i + l)-th assignment cost is at least 

< 

0 , otherwise. 

Number of assignments in the i-th block between the positions 
p + ai + + 1 and p + a, + /3j that are at least . 


With the above notation the formulation of CSA-OMHLP is: 

I G 

{Fuv) min EE EEE “ 1 “ ^Om)^kim ( 1 ^) 

i=l h=2 k£AeeAmeA 

s.t. Constraints : (0) - CU), 

I I 

^a,u^h + ^Vih + ai+i>Y^ Xjk, yh = 2,...,G ( 16 ) 

i=l i=l j£A keA 

Cjk>C(h) 

Uih>Ui-i^h, yi = 2,...,I, h = l,...,G (17) 

Pi-iUih > Vi-i^h, yi = 2,...,I, h = l,...,G (18) 

‘^ih — f ? • • • ? ^ 1, . . . , dr (f^) 

ni/,G{0,l},Ui/,€Zn[0,/3i], Vz = l...,/, h = l,...,G (20) 

Xjfc € {0, l},Sfcrm > 0, yj,k4,m^A. (21) 


The objective function (fT5)l is a reformulation of ([3]) substituting the u variables by the 
new u, V variables and the vector 7 , taking advantage of the A vector properties. Constraints 
(I16p ensure that the number of sites that support a shipping cost to the first hub greater than 
or equal to C(/j) is either equal to the number of allocations with a cost at least C(/j) whenever 
vjh > 0 or less than or equal to 0 : 7+1 otherwise. Constraints (I17p are sorting constraints on 
the variables u similar to constraints m, and constraints (fTSP - ffToP provide upper and lower 
bounds on the variables v depending on the values of variables u. 

The main difference between and F^v is that all Uih variables associated with blocks 
of zero A-values are removed, and those associated with each block of non-null A values 
are replaced by 2 x G variables. Therefore, overall we reduce the number of variables by 
(n - 21) X G. 
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Note that in Formulation Fu, the family of constraints that links covering variables (vari¬ 
ables u) and the allocation variables (variables x), i.e. (1121) . is given with equalities. This fact 
implies that the actual dimension of the feasible region in the space of u and x variables is 
smaller than the one that we were currently working on. This is exploited in the new formula¬ 
tion. Indeed, we reduce the number of variables used in the sorting phase replacing uhy u and 
V. Therefore, the dimension of the feasible region in the space of u, v, x variables has smaller 
dimension. In addition, the constraints that link sorting and design variables, namely (I16p . 
are given as inequalities. This new representation, although valid for the problem, induces 
some loss of information in that it does not allow us to take full control of the exact number 
of allocations at some specific cost. This does not affect the resolution process but influences 
the derivation of valid inequalities. 

Finally, for those cases where /Ij = 1 we observe that Vih = Ui^. This set of constraints 
whenever valid, was added to reinforce the formulation. 

Example 2.1 To illustrate how the F^v versus F^ formulations work, we consider the follow¬ 
ing data. Let A = {1,..., 6} be a set of sites and assume that we are interested in locating 
p = 2 hubs. Let the cost and flow matrices be as follows: 


( 

0 

14 

15 

16 

15 

9 

\ 


( 

0 

15 

2 

8 

11 

2 

\ 


5 

0 

7 

2 

19 

16 




19 

0 

1 

16 

20 

7 



16 

5 

0 

7 

1 

19 


, IV - 


3 

9 

0 

3 

11 

16 



12 

1 

10 

0 

13 

1 




7 

2 

5 

0 

14 

5 



1 

9 

9 

15 

0 

2 




15 

4 

20 

4 

0 

1 


V 

8 

10 

16 

8 

4 

0 



1 

12 

4 

7 

11 

18 

0 

/ 


Therefore, C(.p the sorted vector of c, is in our case 

C(.) = [0,33,42,44,88,126, 208, 210, 294, 315,330, 342,396,416,429,441, 520, 532, 570, 608,660, 
672,798,832,1008,1197]. Hence, G = 26. Let A = (0,1,0,0,1,1), p = 0.7, 5 = 0.9, and the 
capacity constraints vector b = (119,119,113,145,149,140). The optimal solution opens hubs 
4 and 6. The allocation of origin sites to first hub is given by the following values of the 
variables x (see Figure[Jf): 

xie, = X24 = a;34 = X 44 = x^q = xe,e, = 1 - 

Analogously, the allocation of first hubs to final destinations are given by the values of the 
non null variables s. Thus, the flows considering as first hubs 4 and 6 are (see Figure{J\for a 
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graphical representation of the delivery paths): 

S 442 = 11, S 443 = 6, S 444 = 19, S461 = 29, S465 = 45, S466 = 28; 


®642 — 23, Sg44 — 23, Sggi — 27, Sggs — 29, 5665 — 29, 5666 — 3. 



Figure 1: Illustration of Example 12.11 Left figure represents the allocations of sites to their corresponding 
first hubs. Right figure represents the flow pattern to the final destinations from the first hubs: 4 (dashed 
lines) and 6 (grey lines). 

Moreover, the covering variables Uih are given below. Due to their structure, we only report 
for each i the last one and first zero occurrences since they characterize the remaining values. 

t = 1 1-^ till = 1, Ui2 = 0 i = 2 U21 = 1,U22 = 0 t = 3 1-^ Uss = 1, n36 = 0 

t = 4 1 -^ u^Q = 1, U 47 = 0 z = 5 !-)• tisg = 1, n 5 ,io = 0 i = 6 ^ 6,12 = 1, tie ,13 = 0. 

The first two assignments are done at a cost C(i) = 0, corresponding to the two hubs 

(p = 2). The next assignment has been done at a cost C( 5 ) = 88, since U 35 = 1 and = 0, 

and so on. The rest of assignments costs are then C(g) = 126, c^g) = 294 and C( 42 ) = 342. 

Hence, the overall cost of this solution is 
G 

EE EEE {pLCu + Scem)sMm = 636 + 1500.8 = 2136.8. 

i£Ah=2 keAeGAmeA 

In addition, to illustrate how the formulation F^v is related with F^, we also include the 
solution of the covering variables Uih and Vih: 

I = 1 GA UlQ = 1, Ul 7 = 0 

I = li-A Vu = ... = Vig = 2; Pi^io = ... = t;i,i2 = 1; fi,i3 = 0 
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'^i^h — 3 

4 

5 

6 


C(l) 

C(2) 

<^(3) 

C(4) 

= (5) 

= (6) ■ 

• C(24) 

= (25) 

= (26) 
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0 
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0 
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Figure 2: Variables and lambda vector of Example 12.11 


Note that we have only one block of repeated non-null elements of the X-vector, so I = 1. 
(See the right part of FigurelM) The number of zero entries between two blocks is a\ = 2, and 
the number of elements in the 1-st block of non-null elements is /3i = 2. Furthermore, 71 = 1 
is the repeated value in the 1-st block. 

‘=( 1 ) ^( 2 ) *^( 3 ) • • • ‘^( 9 ) ^( 10 ) *^( 11 ) ‘^( 12 ) ^( 13 ) ‘^( 14 ) ‘^( 16 ) ^( 17 ) • • • ‘=( 25 ) ^( 26 ) 


ui,h = (l 1 1 ...1 1 1 1 1 1 1 1 0 ...0 0 ) 

vi,h = (2 2 2 ...2 1 1 1 0 0 0 0 0 ...0 0 ) 

The variable ui^h points out to the row p + ai + 1 = 5 of the original variable Ui^h- Whereas 
the variable vi^h accounts for the number of assignments between the positions p + oi + 1 = 5 
and p + ai + /?! = 6 of Ui^h that are at least C(^h)- (See Figure\^) 

From Ui^h, we know that the 5-th assignment cost is £( 9 ) and the 6 -th assignment cost is 
C(i 2 ) • For this reason vi^h is equal to 2 up to column 9, this is the number of assignment costs 
greater than C(g). Being this number equal to 1 from h = 16 to h = 12, and zero for the 
remaining columns. 

Applying this formulation, the overall reduction in the number of variables is (n —2/) xG = 
104. The rest of variables x and s remain the same, and again the overall cost of this solution 
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is 

I G 

EE 7i(c(/i) - Gh-i))^ih + X] X] X] = 636 + 1500.8 = 2136.8. 

1=1 h=2 kaAlaAm^A 

3 Strengthening the formulation 

3.1 Variable fixing 

Next, we describe some preprocessing procedures that we have applied to reduce further 
the size of formulation F^v We present a number of variable fixing possibilities for the set 
of variables u and v which are useful in the overall solution process. The variable fixing 
procedures developed in this section are based on ideas used in [MIET] and taking advantage 
of the capacity constraints. Indeed, we are adding the reinforced effective capacity constraints, 
as well as some surrogated version of constraints (I24jl since in this aggregated form they give 
better running times. The preprocessing phase developed in this paper also provides new 
upper and lower bounds on the v variables. The percentage of variable reduction obtained by 
these procedures can be found in Tables [U and [2] (column named as ’Fixed’). 

Before describing these procedures for fixing variables, the following simple arguments 
allows us to fix some variables: 

1. First, Cjj = 0 Vj € A, i.e., C(i) = 0. Moreover, every origin where it has been located a 
hub will be allocated to itself as a first hub. 

2 . Second, Cjk 7 ^ 0 if and only if j 7 ^ k, i.e., any non-hub origin is allocated to a first hub 
at a cost of at least C( 2 )- 

Therefore, since in this formulation the first p-allocations are considered only implicitly, we 
can fix un = 1, vn = j3i as well as Ui 2 = 1, Vi 2 = /3j, Vf = 1,..., /. 

3.1.1 Preprocessing Phase 1: Fixing variables to the upper bonnds for the 

formulation with covering variables strengthen with capacity constraints. 

Due to the definition of the variables in formulation F^v , one can expect that Uih = 1 whenever 
i is large and h is small to medium size because this would mean that the {p + ai + -|- 1)- 
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th sorted allocation cost would not have been done at cost less than C(^h)- The reader should 
observe that an analogous strategy applies to the variables v since their interpretation is 
similar, but in this case the values of Vih would be fixed to For the cases where it is not 
possible to fix the corresponding u-variable, it could be still possible to establish some lower 
bounds as we will see later. 

Next, to fix variables Uih and Vih for i = 1,... ,1, h = 1,... ,G, we deal with an auxiliary 
problem that maximizes the number of variables that may assume zero values, satisfying 
Cjk < Gh-i) the capacity constraints. For any h = 1,...,G and j,/c G A such that 
Cjk ^ C(h—1) fof 


jk 


= < 


1 , 


if origin site j is assigned to hub k 


( 22 ) 


j^O, otherwise. 

To avoid possible misunderstanding in the cases where variables are not defined, i.e. when 
Cjk > C(/i-i), we can assume that Zjf. := 0. 

For a given h, we introduce the effective capacity of a hub fe at a cost at most C(/i-i) as. 


:= 


:= mm < 


bk^ 


seA 

Csfc<C(/i_l) 


IT, 


(23) 


Indeed, the capacity of a hub k is always lower than or equal to bk. In addition, if we restrict 
ourselves to the nodes served at a cost of at most then the actual capacity to cover 

this set should be lower than Xlseyt-c k<c(h i) gives us the expression of b^~^. 

The optimal value Pi (h) of the following problem fixes the maximal number of allocations 
that may be feasible at a cost of at most c^k-i)- 


Pi{h) := 


max z!^k 

j,k^A: 

Cjk<C(h-l) 

s.t. Y e 

Y ^jk ^ ^Vk, yk£A 

j€A:Cjk<C(^h-l) 

Yy^-p 

k£A 


( 24 ) 
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ykeA 

j&A:Cjk<C(^h-l) 

z^f,, Vk G {0,1}, Vj, k€A, \/h = l,...,G. 

Then, depending on the value Pi{h) we can fix some variables to their upper bounds. Let 
us denote by ii{h) € {1 ,..., the index such that 


P + < Pi{h) < p + 


Uih = i,Vih = A, * = n{h), 

✓ ' 

Vn{h),h >P + aii(ft) + ^h{h) -Pi{h), 

< > 

Uih = l, Uj/i =/3i, i = + 1,... 

\ / 


if Pi(/i) <p + aii(fe) + Ai(fe)-i 
otherwise. 


3.1.2 Preprocessing Phase 2: Fixing variables to their lower bounds 


Following similar argument to the previous subsection, one can expect that many variables u 
and V in the top-right hand corner of the matrices of variables u and v , respectively, will take 
value 0 in the optimal solution. Indeed, Uih = 0 means that the (p-|-a* -|- /Sj_i + l)-th sorted 
allocation cost is less than which is very likely to be true if h is sufficiently large and i is 
small. Note that an analogous strategy, applies to the variables v since their interpretation is 
similar. For the cases where it is not possible to fix the corresponding u-variable it could be 
still possible to establish some upper bounds. 

For any h = 2,... ,G, j,k £ A such that Cjk > ^{h-i) define variables (122]) • To 

avoid possible misunderstanding in the cases where variables are not defined, i.e. when 
Cjk < we can assume that := 0. Using these variables, the formulation of the 

problem that maximizes the number of non-fixed allocations at a cost at most C(h-i) is: 

P 2 {h) := max ^ z^k 

j,k^A\ 

Cjk>C(h-l) 

s.t. z^k < Vj € A 

kGA:Cj/s>C(h-i) 

Y z'lk<nyk, yj,k£A 

j&A:Cjk>C(^h-l) 


( 25 ) 
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"^yk^p, 

keA 

Vk € {0,1}, Vj, k£A, V/i = 1,..., G. 

Note that the value P 2 {h) implies that there are no feasible solutions of the original 
problem with less than n — P 2 {h) allocations hxed at a cost at most 
Let 1 < i 2 {h) < / be the index such that 

P + a* 2 W-i + ^i 2 {h)-i < n- P 2 {h) <p + ai 2 (^h) + 

Thus, in any feasible solution of the problem we have that: 

Uih = = 0, i = 1,... ,i 2 (/i) - 1, if n - P 2 {h) <p + ai^(h) + 

^ Ui 2 {h),h = 0, <P + ai^ih) + Pi 2 {h) - (n - P 2 {h)), 

> otherwise. 

Uih = 0, Vih = 0, i = 1 ,... ,i2{h) - 1, 

V V / 

Note that whenever n — P 2 {h) = p then there is nothing to fix and therefore no variables are 
set to zero in column h. 

3.2 Valid Inequalities 

In order to strengthen formulation F^v we have studied several families of valid inequalities. 

In fact, taking advantage of previous experience on the non-capacitated version of the problem 
we have borrowed a first family of valid inequalities that are very simple and that have proven 
to be effective in different ordered median problems with covering variables [36^ I37| . This 
family is 



.h = l,.. 

.,G-I, 

(26) 


,.h = I,.. 

.,G-I. 

(27) 


Since, these families are straightforward consequence of the definition of variables u and v 
they have been included in the original formulation. 

In the following, we describe several alternative families of valid inequalities: three sets of 
inequalities, (f28]l . (f29]l - (l3^ . and (I33|l - (l36]) . based on the combination of ordering and capacity 
requirements and two more sets, (1371) and (f38l) - (l3^ . that do not use capacities. 
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3.2.1 First family of valid inequalities: Valid inequalities based on capacity I 

We can add to this model several families of valid inequalities based on capacity issues that 
help in solving the problem by reducing the gap of the linear relaxation and the CPU time to 
explore the branch and bound search tree. 

Observe that the capacity of the set of hubs that may be used to assign origins in ^4 at a 
cost at most is given by 

^k~^^kk 

kGA 

where is the effective capacity at a cost at most defined by (I23h . Recall that, 

although the capacity of a hub k is always lower than or equal to bk, when we restrict to the 
nodes served at a cost of at most then the actual capacity to cover this set should be 

lower than b^~^. Making use of the above observation, we can add the following family of 
constrains as valid inequalities 

WjXjk < bl~^Xkk yh = 2,. .. ,G,k e A (28) 

j'eA 

Cjk<C(h-l) 

which enforces that all the flow sent from origin-hubs at a cost at most cannot exceed 

the effective capacity at that cost. Observe that in the case where b^~^ takes the value 
(11) dominates (I28p . but in the case where b^~^ = X^seA-c k<c(h i) becomes 

J2jGA-.cjk<e(h-i) < 'EjeA:cjk<eiH-i) ^i^kk- Observe that this last valid inequality is an 

alternative surrogation, with capacity coefficients, of constraints Xjk < x^k that although valid 
do not appear in the model because of its large cardinality. This new form of aggregation has 
provided good results in the computational experiments. 

3.2.2 Second family of valid inequalities: Valid inequalities based on capacity II 

This section introduces another family of valid inequalities based on capacity issues that help 
in solving the problem. In order to present these new valid inequalities based on capacity 
requirements we introduce some new notation. Assume without loss of generality that Wi < 
Wj+i for z E A \ {n} and let Wj = Yli=i hUr and S'fc := {z E A : z < fc} for fc E A be a given 
set of origin sites. 
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In case that the effective capacity at a cost at most C(^h-i) is not sufficient to cover the 

n 

demand of 5^, i.e. ~ i® 1®®® than X]s=i ^4(1 “ Xgs), then at most k — 1 

j=i 

origins of Sk can be allocated at a cost lower than or equal to This argument can be 

applied for each h to the corresponding C(h-i) value. Moreover, we have chosen this particular 
structure of consisting of the k origins (nodes) with the /c-smallest flows, because given 
a fixed amount of flow, the maximal cardinality set of origins, such that the overall flow 
originated in this set is lower than or equal to this amount, is provided by a set for some 
k £ A. Therefore, since we are dealing with the worst cases, it allows us to fix some variables u 
and V through the following valid inequalities. For each h = 2,... ,G we obtain the following: 

• If ai-i + < k <ai + lSi_i + I, for some i G {I,..., /}, namely if the index of the 

last element, /c, that defines Sk lies in the i-th block of null elements in the A vector, 
then 

k 

WkUih + (b)-^ - Wj) Xjj >Y,Ws{l- Xss). (29) 

jeA s=i 

Observe that the above inequality amounts to a disjunctive condition: either the effective 
capacity at a cost at most C(/i-i) is enough to cover the demand of the k smallest flows 
from origin sites or the i-th sorted cost allocation must be assigned at a cost at least 
C{h)- 

• If cfi + /3j_i + I < A: < Oj + /3j, and namely if the index of the last element, k, that 
defines Sk lies in the f-th block of non-null elements, then 

k 

WkVih + Y, (hY - > Y Ws{l - Xss). (30) 

jeA s=i 

In this case, the inequality is similar to the previous one but written in terms of the 
variables v that allow to control the capacity whenever k falls within a block of non-null 
elements in the A vector. 

Remark 3.1 Recall that if k > Pi{h) — p, variables Uih and Vih have been already fixed by 
the Preprocessing Phase 1, and for the above inequality to be effective k < Pi{h) — p, or 
equivalently, i < ifih) . 
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Based on the same arguments we can add a larger family of valid inequalities built on 
arbitrary sets of origin sites. Let 5 be a set of origin sites, and suppose that ^4^ = Xlses 
satisfies Wk < ^5 < Wk+i , 

• If oii-i + j3i_i < k <ai + /3i_i + 1, for some z € {1,..., /}, and k < Pi{h) — p then 

k 

Asu,h + ^ ^ - xss). (31) 

j£A s=l 

• If aj + /3j_i + 1 < k <ai + Pi, and k < Pi{h) — p then 

k 

Asvih + ^ (h)-^ - VL,(1 - Xss). (32) 

j&A s=l 

Now, assuming a more general case and for an improvement of the above inequalities (|29l) 
and (l30)) . for a given /i € {2,..., G}, and k < Pi{h) — p, k {1 ,... , a/+i + Pi)}, let 


Ms := min c,*. 
i(^s)GA 

Ms is the minimum allocation cost to s as an open hub. In other words, no allocation to hub 
s is possible at a cost less than Ms, except in case s were a hub itself. We shall call this value 
the empty radius of s. 

Define s{h — 1, k) to be the index of the sorted sequence of elements Ws such that there 
are exactly k elements Ws with s < s{h — 1, k) and Mg < namely s{h — 1, k), is the 

index such that 

|{s : s <s{h-l, k); Mg < C(/i_i)}| = k. 

Then it holds that. 


• If ai-i + Pi_i < k < Oi + Pi— I + 1, for some i € {1,..., /}, and k < Pi{h) -p 


s{h—l,k) 

w s^n-i,k)-Y.^^ 

s = l 


Uih + (j>j — 

j&A 


s(/i—1,/c) s(/i—1,/c) 

^ ^ “ ^ss) ~ ^ ^ 1,/c) “ ^s)^ss 

S=1 S=1 

^s<C(h-l) Afs>C(/i_i) 


( 33 ) 
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The above inequality is also a disjunctive condition that “reinforces” the family of valid 
inequalities (I29p . It states that if the effective capacity at a cost at most C(h-i) is 
not enough to cover the flow sent from origin sites that are not hubs and that can be 
allocated at some costs less than or equal to C(^h-i) then some of the origin sites with 
allocation costs less than C(/i-i) must be assigned at a cost at least C(/jp We observe 
that the use of u variables in the left-hand side of the inequality is due to the fact that 
k falls within a block of null elements in the A vector. A similar inequality also holds 
when k falls within a block of non-null elements as shown below. 


If ai + I3i_i + I<k<ai + j3^, and k < Pi{h) — p then 


/ 


W 


s(h—l,k) 




jGA 


s(h—l,k) 

S=1 

s(h—l^k) s{h—l,k) 

^ W,(l - Xss) 


(34) 


S=1 S=1 

^s<C(h-l) Ms>C(^h-l) 

This inequality is similar to the previous one whenever the index k falls within a block 
of non-null elements in the A vector. 


Finally, as the index s{h — I,k) should be greater than or equal to k, we can split the 
above equations, (1331) and (1341) . into s{h — l,k) — k equivalent inequalities. This is, for any 
t = k,..., s{h — 1, k), define s{h — 1, k,t) to be the index of the sorted sequence of elements 
Wg such that 


|{s ■, s < s{h — 1, k,t)', Mg < C(/i_i)}| + t — s{h — l,k,t) = k. 

Then it holds that. 


• If ai-i + <k < ai + -|- 1, for some i € {1,...,/}, and k < Pi{h) — p 


/ 

Wt 

V 


- t 

s=l 

Ms>c^^_i‘^a.nd s<s(h—l,k,t) j 


Uih + (bj — 

j£A 


t t 

Y, ^ {Wt-Wg)xgg 


s = l 

Ms<C(^_i)or s>s{h—l,k^t) 


S = 1 

Msand s<s(h—l^k^t) 


( 35 ) 
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• If ctj + /3j_i + 1 < k < ai + Pi and k < Pi{h) — p, then 

(_ 

w, - Y, 

S = 1 

t t 

^ Ws{l-Xss) - Y. iWt-Ws)Xss. (36) 

S=1 S=1 

Ms<C(^h—i)^^ s>s(h—l,k,t) Ms>C(^iT^_i),s<s{h—l,k,t) 

Observe that if t = s{h — 1, k) then s{h — 1, k, t) = s{h — \,k). Thus, the families of valid 
inequalities ()35p and (I36p include as particular instances the families (|33l) and (|34p . 

3.2.3 Third family of valid inequalities: Disjunctive implications 

The third family of valid inequalities, directly borrowed from m, state disjunctive implica¬ 
tions on the origin-first hub allocation costs. They ensure that either origin site j is allocated 
to a first hub at a cost of at least C(/j) or there is an open hub k such that Cjk < C{h)- This 
argument can be formulated through the following family of valid inequalities: 

Z] + z Xkk>l, Vj G A = (37) 

k&A:Cjk>C(h) kGA:Cjk<C(^h) 

3.2.4 Fourth family of valid inequalities 

Using the definition of the variables u and v, we establish a lower and an upper bound of the 
number of feasible allocations at a cost Observe that using the family of constraints 

m for the original formulation, the exact number of allocations done at a cost C(^h-i) is 
given by ~ Ui,h)- However, since in formulation F^v the number of variables 

has been considerably reduced, some information is lost. In particular, we cannot keep under 
control with this new formulation the exact number of allocations at a cost Indeed, the 

counterpart to equalities (I12p in formulation F^v is the family of constraints (I16p . Therefore, 
we are only able to give a lower and upper bound of this number of allocations. These lower 


Vih + (y^j — 

jeA 
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and upper bounds are formulated, respectively, by the following two families of constraints: 

I I 

P ^ 1 '^ih) “1“ ^ ^ 1 Uih^, V /l — 2, . . . , G, (38) 

i=l i=2 

I I 

< - Vi^h) + ^ - Ui^h) + (1 - Uih)ai+i + 

i=l i=l 

I-l 

“1“ ^ ^ Oii+l ^ih)j h 2, . . . , G. (39) 

i=l 

The first sum in the right hand side of both families gives the exact number of allocations 
at cost C(h-i) ill ths positions corresponding to non-null blocks of the vector A. However, the 
second sum in the right hand side of constraints (13811 provides a lower bound of the number of 
allocations at cost Cf^h-i) in tii® positions corresponding to the null blocks. In the same way, 
the 2nd to the 4-th sums in (|39p provide an upper bound on the number of these allocations. 

4 Computational Results 

The formulations given to the CSA-OMHLP with the corresponding strengthening and pre- 
proccesing phases, described in this paper, were implemented in the commercial solver XPRESS 
IVE 1.23.02.64 running on a Intel(R) Core(TM) i5-3450 CPU @3.10GHz 6GB RAM. 

The cut generation option of XPRESS was disabled in order to compare the relative 
performance of the formulations cleanly. 

Eor this purpose we use the AP data set publicly available at http://www.cmis.csiro.au/or 
/hubLocation (see US]). As in previous papers on the field related to the uncapacitated version 
of this problem, we tested the formulations on a testbed of five instances for each combination 
of costs matrices varying: (i) n in {15, 20, 25, 28, 30} (ii) three different values of p depending 
on the case and (iii) p = 0.7, 5 = 0.9/i and six different A-vectors. These A-vectors are the well- 
known Median A = (!,...,!), Anti-(A:i -|- /c 2 )-trimmed-mean A = (1, 1,0,..., 0,1, 1), 

{ki + /c 2 )-Trimmed-mean A = (0, .^.7,0,1,... , 1,0, with ki = k 2 = [0.2n], Center A = 

(0,... ,0,1), and A:-Centrum A = (0,... ,0,1, 1) with k = [0.2n]. As well as a {0, l}-blocks 

A-vector (three alternate {0 — l}-blocks of lambda weights, i.e. A = (0,... , 0,1,..., 1,0,... , 0, 


Z] Z 

j£A _ keA 

Cjk=C(^h-l) 

Z Z 

jeA keA 
Cjk=Lh-l) 
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□ F, ^ F„„ □ F„„+Pi2 ^ F„,+ Pi2+VI(27),(34)-(35) H F„„+Pi2+VI(27),(34)-(38) 




Figure 3: Summary of computational results carried out in the paper. 

1,..., 1,0,... , 0,1,... , 1)). Therefore, for the each combination of n, p and A we have tested 
five instances. This is, a total number of 450 problems have been used to test the performance 
of the proposed models. 

The capacities were randomly generated in [minPFj, (1/2) > Wi]. This generation proce- 

i&A 

dure does not ensure in all cases feasible instances, as capacity constraints can be very tight 
for problems with a low number of hubs (n = 3,5). Overall, in our experiments we got, ini¬ 
tially 10 infeasible instances out of 75 (13.3%). These instances were replaced by new feasible 
ones (generated with the same capacity structure). Hence, the reader may observe that the 
generation procedure gives tight capacity constraints. 

First of all, and for the sake of readability, we present in Figure [3] a summary of our 
computational results. A full description of those results is also included in Tables [T](6j 

Figure [3] shows average results for each one of the considered problem types (different 
A-vectors). The left chart refers to average CPU times and the right chart to the number of 
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explored nodes of the B&B tree. Both charts contain the same number of blocks standing for 
the different types of A-vectors plus and additional block (Total Avg.) for the consolidated 
average of all A-vectors. Each block compares the behavior of the different formulations 
{Fu,Fuv) and their strengthening (variable fixing and valid inequalities). The heading explains 
the meaning of bars within each block: and Fu^ stand for the corresponding formulations; 

FuvTPi 2 when the two preprocesses Pi{h) and P 2 {h) are applied; and E„,;+Pi 2 +Vl when the 
Valid Inequalities are added as well, denoted with their corresponding references. 

Analyzing this figure we observe the improvement obtained with F^v and its strengthening 
as compared with Fu or even F^v alone. Actually, the overall reduction in running time with 
respect to the initial formulation is above 89%. 

Next, focusing in the best model, namely K,„ + Pi9+VI (l28p . (l35p - (l39p . we observe that the 
most time consuming problem is the one with A-vector given by {0, l}-blocks with a significant 
difference with respect to the remaining A-vectors. The second most time consuming problem 
corresponds with the Antitrimmean. Similar conclusions, regarding the number of nodes, are 
obtained looking at the right chart of Figure [3l It is worth mentioning that although the 
formulation F^ provides the worst computational times, it is the one that reports the lowest 
number of nodes in the B&B tree. This fact is explained because F^ has a larger number of 
variables than which results in more difficult LP relaxations in each node of the B&B 
tree. 

In spite of that Figure [3] shows the general overview of our computational results, we also 
include in the following a more detailed analysis based on Tables (H-El 

Tables [1] and [2] report the results of the formulations and different preprocessing phases 
developed in this paper for the CSA-OMHLP. The first column of these tables describes the 
different A-vectors in the study, the second and third columns report the size of the instances 
and the number of hubs to be located, respectively. The following three columns correspond to 
some of the computational results obtained by solving the CSA-OMHLP with F^ formulation. 
The next three columns correspond to Formulation P„^, and the rest of the columns to the 
latter formulation plus the two preprocessing procedures, i.e Ptjt;+Pi+P 2 - Columns RGAP, 
Nodes and Time stand for the averages of: the gap in the root node, number of nodes in the 
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B&B tree and the CPU time in seconds; the time was limited to two hours of CPU. To obtain 
a general idea of the comparisons among these averaged values, for the results in the column 
Time and for different formulations and/or valid inequalities applied, we have accounted the 
value 7200 seconds for those instances that exceed the time limit. A superindex in their 
corresponding averaged time value states the number of instances exceeding the CPU time 
limit; in the same way, the values used to computed the average of the column Nodes have 
been the number of nodes of the B&B tree when the CPU time limit was reached. 

The column Fixed gives the percentage of variables that have been fixed after the Prepro¬ 
cessing Phases 1 and 2. Column Cuts provides the number of the lower and upper bounds 
over the variables v added to the model, after running the corresponding preprocessing phases. 
Finally, T. prep, reports the CPU time in seconds of the corresponding preprocessing phases 
and column T. total reports the overall CPU time in seconds to solve the problem includ¬ 
ing the corresponding preprocessing phase. The rows Average provide the averaged results 
among all the tested instances for each problem and TOTAL with respect to the overall set 
of instances. 

Tables [T] and [2] show, as a general trend, that the CPU time increases similarly, for all 
choices of the A-vector, with the size of the instances. We can see that 69 instances required 
more than two hours to be solved using Formulation however the remaining two analyzed 
ways to solve this problem were able to solve all the instances within the CPU time limit, 
except in a few cases for the Antitrimmean and {0, l}-blocks problem types. Regarding the 
running times for the different types of problems, we can see that Anti-TrimMean and Blocks 
have been the problems that need more time to be solved. In any case, we observed that 
there is a considerable reduction of the running times from the original formulation to the 
improved formulation, and after applying the preprocessing phases. From tables [T] and [2l 
one can remark that, Formulation F^v with Preprocessing Phases 1 and 2 provides better 
results, with a reduction of around a 86% of the time with respect to Formulation F^ (taking 
into account that the latter formulation was not able to solve all the studied instances before 
the time limit was exceeded). Moreover, the reduction of the running times with respect to 
Formulation F^v (without preprocessing phases) is around 32 %. 


5 CONCLUDING REMARKS 


27 


As for the comparison between RGAP’s, we observe that the average gap of the linear 
relaxation after preprocessing reduces around 32% from the original formulation (Formulation 
Fu) to Formulation Fu^ with Preprocessing Phases 1 and 2. In any case, it is also worth 
noting that the gap from the original formulation to the improved formulation (without any 
preprocessing phase) is reduced around 9%, what implies that even though this improved 
formulation uses much less number of variables, it provides a better RGAP. 

Tables [3] and m present several improvements to Formulation with Preprocessing Phases 
1 and 2. In particular, the second, third and fourth blocks of columns summarize the results 
when the combination of valid inequalities ([28]) - (l30]) : ([28]) . (l33]) - (l34l) and (I28]) . (p5]) - (l36|) are 
added to this formulation, respectively. In general, we can see that the latter combination 
provides the best results. In particular, there is an improvement of 7 % of the running times 
with respect to Formulation F^y with Preprocessing Phases 1 and 2. Moreover, with respect 
to the number of nodes the improvement is around 15%, but the RGAP is similar for all the 
analyzed reinforcements. 

Since the best behavior observed was obtained with reinforcement given by (1281) . (I35l) - (f36]) . 
in the rest of our tests we have used this configuration to make further strengthening. Tables 
[5] and [6] present several improvements to the Formulation Fyy with Preprocessing Phases 1 
and 2 and valid inequalities (I28p . (p5|) - (l36]) . In particular, the second, third and fourth blocks 
of columns report the results after including the family of valid inequalities (I37p . (I37p - (l38]) 
and (j37p - ([M]) . respectively. These tables show that the best results are obtained when all 
valid inequalities, (I37l) - (l39p . are added to the current configuration; with an improvement of 
18% in the running time, 28% in the number of nodes and 12 % in the RGAP. 

The overall conclusion of our experiment is that in order to solve GSA-OMHLP the best 
combination of formulation and strengthening is to use K,„+i-i9+ (|^ + (l35P + (l36l) + (l371) + (l39p . 
This configuration allows to solve medium size instances within 10 minutes of CPU time. 

5 Concluding remarks 

This paper can be considered as an initial attempt to address the capacitated single-allocation 
ordered hub location problem. The formulations, strengthening and preprocessing phases de- 
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veloped in this paper provide a promising approach to solve the above mentioned problem 
although so far only medium size problems are reasonably well-solved. Thus, this work 
opens interesting possibilities to study-and-develop ad-hoc solution procedures that allow us 
to consider larger size instances of this problem. Moreover, it also points out the possibility of 
developing heuristics approaches that will give good solutions in competitive running times. 
All in all, this paper shows the usefulness of using covering formulations and their correspond¬ 
ing strengthening for solving capacitated versions of ordered hub location problems. 
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Table 1: Computational results for Formulations Fn vs F^v with preprocessing. 
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Table 2: Computational results for Formulations Fn vs F^v with preprocessing(II). 
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76 


20 

3 

21,09 

606 

113,0 

20,86 

3721 

30,7 

17,58 

5777 

24,7 

3,2 

16,1 

219 


20 

8 

8,55 

402 

67,3 

8,55 

12670 

41,3 

7,46 

8064 

20,4 

2,2 

4,8 

116 


20 

10 

4,70 

97 

33,2 

4,32 

3590 

8,4 

4,00 

4585 

9,0 

1,7 

3,7 

122 

Vi 

25 

3 

21,51 

4850 

983,7 

21,27 

5637 

114,4 

16,45 

4087 

59,1 

8,9 

27,5 

484 

u 

o 

u 

CQ 

25 

8 

9,99 

1675 

429,2 

9,94 

13400 

115,0 

8,12 

12050 

89,3 

5,2 

17,4 

260 

25 

10 

17,57 

6874 

1021,0 

17,40 

46270 

216,9 

10,01 

40870 

177,5 

5,2 

7,4 

268 


28 

3 

22,45 

5627 

2887,0 

22,45 

9111 

472,4 

17,16 

8569 

267,1 

13,4 

48,3 

482 


28 

8 

10,44 

14490 

2687,0 

10,25 

96840 

1515,0 

9,32 

54910 

745,0 

7,7 

24,1 

357 


28 

10 

9,39 

13040 

3181,0^^^ 

9,26 

233200 

2683,0 

8,67 

184700 

1601,0 

6,7 

15,4 

218 


30 

3 

33,18 

10180 

5755,0^2> 

33,04 

10250 

753,4 

22,39 

12840 

703,5 

18,9 

56,4 

664 


30 

8 

21,41 

15590 

6438,0^^^ 

21,46 

201400 

6587,0 

12,95 

193300 

5016,0^^^ 

12,9 

24,5 

511 


30 

10 

11,90 

18420 

4144,0^^^ 

11,66 

283700 

5113,0 

9,45 

256600 

3402,0^^^ 

11,51 

23,28 

465 


Average 

15,80 

6132 

1851,6 

15,60 

61612 

1177,7 

11,60 

52644 

808,1 

6,8 

19,2 

297 

TOTAL 

17,40 

4268 

1966,3 

15,90 

34118 

420,1 

11,80 

26543 

286,5 

6,7 

8,8 

279 
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Table 3: Computational results for FormulationF^^ with Valid Inequalities (I) 





Pi 

+ P 2 - P 

12 

Pl 2 

+V.l.i28t-i30t 


Pi-)+V.l.i2^.{|3^-{|36t 




RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 


15 

3 

15,16 

195 

1,4 

14,21 

199 

1,4 

14,21 

262 

1,5 

14,21 

239 

1,4 


15 

5 

6,78 

177 

0,7 

6,50 

154 

0,7 

6,50 

232 

0,8 

6,50 

230 

0,8 


15 

8 

4,39 

185 

0,4 

4,33 

110 

0,4 

4,33 

157 

0,4 

4,34 

162 

0,5 


20 

3 

16,31 

472 

5,8 

15,66 

379 

7,1 

15,66 

540 

7,4 

15,66 

560 

7,3 


20 

8 

5,49 

1122 

3,8 

4,98 

981 

4,1 

4,98 

1189 

4,4 

4,98 

936 

3,9 


20 

10 

2,59 

560 

1,6 

2,33 

658 

1,9 

2,33 

654 

1,8 

2,33 

857 

2,4 


25 

3 

15,21 

408 

18,5 

14,79 

502 

21,3 

14,79 

490 

20,9 

14,79 

390 

19,9 

1 

25 

8 

5,95 

6177 

33,6 

5,39 

1772 

16,8 

5,39 

2994 

21,6 

5,39 

3301 

26,5 


25 

10 

6,31 

4686 

26,0 

5,97 

2748 

16,1 

5,97 

3141 

18,5 

5,97 

2734 

16,8 


28 

3 

16,01 

2246 

76,6 

15,72 

1968 

72,6 

15,72 

1871 

71,5 

15,72 

2055 

71,9 


28 

8 

8,15 

10560 

97,4 

7,53 

5937 

62,5 

7,53 

5999 

66,2 

7,53 

6398 

70,4 


28 

10 

6,78 

32680 

243,1 

6,48 

21080 

155,2 

6,48 

16800 

136,4 

6,48 

17930 

125,4 


30 

3 

20,92 

2564 

127,8 

20,60 

2300 

132,9 

20,60 

2797 

149,9 

20,60 

2770 

142,3 


30 

8 

10,86 

36790 

437,3 

10,35 

31990 

448,0 

10,35 

25920 

338,3 

10,35 

25060 

332,7 


30 

10 

6,98 

20130 

207,6 

6,75 

17230 

167,8 

6,75 

15720 

160,4 

6,75 

18130 

179,8 


Average 

9,90 

7930 

85,4 

9,40 

5867 

73,9 

9,40 

5251 

66,7 

9,40 

5450 

66,8 


15 

3 

17,62 

645 

2,0 

17,18 

559 

1,9 

17,18 

709 

2,1 

17,18 

556 

2,0 


15 

5 

8,38 

803 

1,4 

8,36 

823 

1,6 

8,36 

615 

1,4 

8,36 

705 

1,5 


15 

8 

6,82 

724 

0,9 

6,82 

808 

0,9 

6,82 

723 

0,9 

6,82 

736 

1,0 


20 

3 

17,44 

1150 

9,2 

17,00 

1471 

11,4 

17,00 

1455 

11,7 

17,00 

1370 

10,9 


20 

8 

7,66 

3009 

7,8 

7,54 

2684 

7,1 

7,54 

2716 

7,6 

7,54 

2709 

7,7 


20 

10 

4,28 

3394 

6,4 

4,28 

2663 

6,1 

4,28 

2192 

5,1 

4,28 

2574 

5,7 

1 

25 

3 

17,29 

1218 

29,0 

16,95 

1294 

31,7 

16,95 

1314 

33,9 

16,95 

1167 

28,5 


25 

8 

7,91 

6841 

33,31 

7,74 

5315 

30,7 

7,74 

3127 

22,7 

7,74 

5038 

30,4 


25 

10 

9,81 

25390 

102,1 

9,76 

11640 

55,2 

9,76 

12560 

61,8 

9,76 

13830 

66,2 


28 

3 

17,36 

7432 

168,4 

17,10 

6343 

155,4 

17,10 

7032 

164,6 

17,10 

7531 

168,8 


28 

8 

10,67 

34800 

271,1 

10,63 

23080 

204,8 

10,63 

21210 

191,3 

10,63 

20850 

188,2 


28 

10 

9,141 

131900 

832,1 

9,11 

92660 

659,0 

9,11 

81600 

634,5 

9,11 

93700 

705,6 


30 

3 

21,74 

7220 

271,6 

21,47 

3502 

164,5 

21,47 

4853 

204,9 

21,47 

4873 

190,1 


30 

8 

14,48 

273600 

2896,0 

14,36 

183100 

2022,0 

14,36 

204500 

2250,0 

14,36 

258800 

3026,0 


30 

10 

10,56 

109700 

863,2 

10,50 

99160 

837,6 

10,50 

94870 

784,8 

10,50 

103700 

898,1 


Average 

12,10 

40522 

366,3 

11,90 

29007 

279,3 

11,90 

29298 

291,8 

11,90 

34543 

355,4 


15 

3 

16,44 

729 

1,4 

16,19 

677 

1,4 

16,19 

733 

1,4 

16,19 

789 

1,5 


15 

5 

8,58 

646 

0,7 

8,39 

566 

0,7 

8,39 

550 

0,8 

8,39 

489 

0,7 


15 

8 

7,88 

322 

0,3 

7,85 

353 

0,3 

7,85 

349 

0,3 

7,85 

366 

0,3 


20 

3 

17,65 

1022 

6,2 

17,60 

871 

6,5 

17,60 

977 

6,2 

17,60 

929 

6,8 


20 

8 

9,08 

4512 

7,9 

9,00 

4519 

8,3 

9,00 

4508 

7,9 

9,00 

4476 

8,4 


20 

10 

5,97 

1620 

2,7 

5,95 

1402 

2,5 

5,95 

1353 

2,5 

5,95 

1444 

2,6 

S 

25 

3 

16,87 

2081 

25,1 

16,83 

1844 

24,1 

16,83 

1879 

26,8 

16,83 

1486 

22,5 


25 

8 

8,63 

21460 

74,2 

8,40 

13420 

54,4 

8,40 

14300 

56,4 

8,40 

15390 

59,5 


25 

10 

12,28 

1760 

6,5 

12,14 

1106 

4,8 

12,14 

1384 

5,7 

12,14 

1251 

5,6 


28 

3 

17,49 

3800 

69,6 

17,48 

4117 

80,5 

17,48 

4005 

71,0 

17,48 

4458 

80,8 


28 

8 

10,24 

20290 

130,9 

9,97 

15890 

111,4 

9,97 

17340 

119,1 

9,97 

16910 

119,1 


28 

10 

9,17 

14770 

69,3 

8,92 

8775 

41,6 

8,92 

9241 

44,4 

8,92 

10490 

49,7 


30 

3 

22,08 

3302 

124,4 

22,05 

4052 

157,8 

22,05 

4356 

157,0 

22,05 

3925 

147,3 


30 

8 

15,09 

54920 

483,1 

14,99 

45860 

450,8 

14,99 

50140 

486,9 

14,99 

52040 

514,7 


30 

10 

9,36 

107100 

757,0 

9,17 

42140 

316,1 

9,17 

36570 

275,1 

9,17 

36480 

263,3 


Average 

12,50 

15889 

117,3 

12,30 

9706 

84,1 

12,30 

9846 

84,1 

12,30 

10062 

85,5 
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Table 4: Computational results for FormulationF^j; with Valid Inequalities (I) 



n 



P1 + P2 ^ 

P12 

P12 

+V.l.i28t-i30t 


Pi-)+V.I.{l2^.i35t-{|3^ 



RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 


15 

3 

16,90 

489 

0,9 

16,90 

429 

0,9 

16,90 

468 

0,9 

16,90 

701 

1,1 


15 

5 

8,48 

809 

0,9 

8,47 

788 

0,9 

8,47 

800 

1,0 

8,47 

759 

1,0 


15 

8 

8,83 

501 

0,6 

8,83 

495 

0,6 

8,83 

554 

0,6 

8,83 

556 

0,6 


20 

3 

16,72 

1017 

4,8 

16,70 

840 

4,6 

16,70 

1028 

4,7 

16,70 

876 

4,5 


20 

8 

11,58 

3712 

6,2 

11,58 

3185 

5,9 

11,58 

3778 

6,6 

11,58 

4408 

7,9 


20 

10 

8,82 

3019 

5,3 

8,82 

2430 

4,9 

8,82 

2377 

4,9 

8,82 

2044 

4,7 

1 

25 

3 

14,57 

938 

12,3 

14,56 

956 

14,3 

14,56 

821 

12,5 

14,56 

983 

13,6 

1 

25 

8 

11,33 

26520 

87,4 

11,33 

25490 

89,2 

11,33 

26310 

95,2 

11,33 

30540 

107,2 


25 

10 

16,16 

5134 

17,7 

16,16 

3933 

15,3 

16,16 

4357 

15,9 

16,16 

5271 

19,2 


28 

3 

17,56 

4172 

66,8 

17,55 

5275 

75,3 

17,55 

4087 

69,8 

17,55 

5190 

86,6 


28 

8 

15,77 

29450 

170,1 

15,77 

25760 

157,3 

15,77 

30410 

185,7 

15,77 

25740 

160,5 


28 

10 

14,53 

48120 

218,4 

14,53 

49170 

233,4 

14,53 

62630 

298,0 

14,53 

53050 

269,0 


30 

3 

19,75 

6536 

147,3 

19,73 

4810 

124,0 

19,73 

6142 

151,6 

19,73 

5994 

151,6 


30 

8 

14,63 

103700 

786,6 

14,63 

86370 

680,3 

14,63 

113900 

904,5 

14,63 

108300 

965,3 


30 

10 

11,77 

96140 

593,6 

11,77 

128200 

882,8 

11,77 

108700 

766,9 

11,77 

94490 

675,7 


Average 

13,8 

22017 

141,3 

13,8 

22542 

152,7 

13,80 

24424 

167,9 

13,80 

22594 

164,6 


15 

3 

17,60 

588 

1,8 

17,58 

659 

1,9 

17,58 

707 

2,0 

17,58 

783 

1,9 


15 

5 

7,91 

498 

1,0 

7,85 

653 

1,3 

7,85 

547 

1,2 

7,85 

635 

1,2 


15 

8 

4,40 

185 

0,4 

4,33 

110 

0,4 

4,33 

157 

0,4 

4,34 

162 

0,4 


20 

3 

18,23 

1.528 

9,1 

18,22 

1590 

10,6 

18,22 

1287 

10,0 

18,22 

1777 

10,7 


20 

8 

7,12 

2.902 

6,2 

7,09 

3290 

7,8 

7,09 

2568 

6,2 

7,09 

2821 

6,7 


20 

10 

4,04 

3705 

5,9 

4,03 

2973 

5,2 

4,03 

2769 

4,8 

4,03 

3045 

5,4 

I 

25 

3 

17,48 

2693 

32,3 

17,47 

1764 

28,4 

17,47 

1537 

26,6 

17,47 

1639 

28,0 

u 

u 

25 

8 

7,33 

3910 

23,2 

7,30 

3239 

21,0 

7,30 

3488 

24,3 

7,30 

3264 

23,5 

« 

25 

10 

7,34 

12820 

52,6 

7,28 

12730 

51,9 

7,28 

11770 

49,6 

7,28 

11330 

48,8 


28 

3 

17,17 

4971 

111,3 

17,16 

4126 

103,2 

17,16 

4417 

106,4 

17,16 

4146 

108,7 


28 

8 

9,44 

21180 

178,9 

9,39 

17390 

158,4 

9,39 

14370 

140,5 

9,39 

14120 

139,8 


28 

10 

7,63 

35460 

228,2 

7,54 

39050 

276,9 

7,54 

38760 

268,8 

7,54 

34770 

247,2 


30 

3 

21,74 

4588 

187,6 

21,72 

4740 

188,6 

21,72 

4069 

169,9 

21,72 

4437 

194,9 


30 

8 

12,20 

171500 

1830,0 

12,14 

149900 

1737,0 

12,15 

218900 

2710,0 

12,14 

114200 

1343,0 


30 

10 

7,63 

37330 

340,0 

7,54 

41100 

392,6 

7,54 

34690 

332,8 

7,54 

30780 

305,4 


Average 

11,20 

20257 

200,6 

11,10 

18888 

199,0 

11,10 

22669 

256,9 

11,10 

15194 

164,4 


15 

3 

15,97 

1729 

3,8 

15,86 

2184 

4,2 

15,86 

1481 

3,7 

15,86 

2248 

4,4 


15 

5 

8,08 

860 

2,4 

7,98 

1085 

2,8 

7,98 

877 

2,4 

7,98 

877 

2,7 


15 

8 

6,25 

715 

1,1 

6,23 

748 

1,2 

6,23 

772 

1,1 

6,23 

726 

1,1 


20 

3 

17,58 

5777 

24,7 

17,41 

3515 

19,7 

17,41 

4382 

19,3 

17,41 

4007 

21,0 


20 

8 

7,46 

8064 

20,4 

7,40 

5603 

15,8 

7,40 

6327 

16,1 

7,40 

5217 

15,8 


20 

10 

4,00 

4585 

9,0 

3,98 

4308 

8,6 

3,98 

4091 

8,5 

3,98 

3505 

7,8 

CO 

25 

3 

16,45 

4087 

59,1 

16,28 

4489 

67,7 

16,28 

4821 

65,1 

16,28 

4434 

63,9 

CJ 

o 

nJ 

oa 

25 

8 

8,12 

12050 

89,3 

8,07 

8427 

64,9 

8,07 

7338 

54,5 

8,07 

7836 

64,9 

25 

10 

10,01 

40870 

177,5 

9,99 

49940 

199,6 

9,99 

36030 

155,2 

9,99 

35520 

158,8 


28 

3 

17,16 

8569 

267,1 

17,11 

11070 

312,0 

17,11 

5840 

220,4 

17,11 

5934 

210,7 


28 

8 

9,32 

54910 

745,0 

9,14 

47500 

723,0 

9,14 

54830 

749,3 

9,14 

43450 

580,5 


28 

10 

8,67 

184700 

1601,0 

8,61 

124900 

1118,0 

8,61 

149500 

1385,0 

8,61 

140900 

1323,0 


30 

3 

22,39 

12840 

703,5 

22,20 

10580 

572,6 

22,20 

10520 

542,8 

22,20 

9661 

517,0 


30 

8 

12,95 

193300 

5016 , 

12,76 

261900 

6298 , 

12,75 

197700 

4858 , 

12,74 

187400 

4854 , 


30 

10 

9,45 

256600 

3402 , 

9,29 

248300 

3476 , 0 ^^' 

9,29 

215800 

3277 , 

9,30 

253000 

3586 , 


Average 

11,60 

52644 

808,1 

11,50 

52303 

858,9 

11,50 

46687 

757,2 

11,50 

46981 

760,8 

TOTAL 

11,80 

26543 

286,5 

11,70 

23052 

274,7 

11,70 

23029 

270,8 

11,70 

22470 

266,2 
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Table 5: Computational results for FormulationT^,; with Valid Inequalities (II) 





Pl2+V.I.{28),l3S)-l36) 


PiQ4-V.l.i28t.{|3^-{|38t 

PT)+V.l.i28t.{|35t-{|3^ 




RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 


15 

3 

14,21 

239 

1,4 

7,71 

168 

2,3 

7,71 

197 

2,5 

7,63 

210 

2,3 


15 

5 

6,50 

230 

0,8 

4,58 

122 

1,4 

4,58 

117 

1,4 

4,60 

132 

1,6 


15 

8 

4,33 

162 

0,5 

4,04 

166 

0,7 

4,04 

150 

0,8 

4,04 

153 

0,9 


20 

3 

15,66 

560 

7,3 

7,08 

220 

11,9 

7,08 

345 

12,3 

7,16 

281 

12,7 


20 

8 

4,98 

936 

3,9 

4,38 

262 

4,4 

4,38 

600 

5,4 

4,36 

318 

4,3 


20 

10 

2,33 

857 

2,4 

2,24 

354 

3,2 

2,21 

328 

2,9 

2,21 

273 

2,9 


25 

3 

14,79 

390 

19,9 

7,52 

429 

44,5 

7,512 

319 

42,9 

7,52 

393 

42,7 

1 

25 

8 

5,39 

3301 

26,5 

4,03 

521 

17,8 

4,02 

486 

18,5 

4,02 

405 

16,3 


25 

10 

5,97 

2734 

16,8 

5,48 

1177 

15,9 

5,48 

1265 

16,4 

5,48 

1175 

17,3 


28 

3 

15,72 

2055 

71,9 

6,75 

320 

88,1 

6,75 

414 

82,8 

6,76 

322 

85,7 


28 

8 

7,53 

6398 

70,4 

5,15 

1627 

47,6 

5,13 

1294 

42,9 

5,15 

1781 

49,3 


28 

10 

6,48 

17930 

125,4 

4,73 

4394 

57,9 

4,73 

4914 

62,8 

4,70 

3674 

50,2 


30 

3 

20,60 

2770 

142,3 

13,90 

868 

143,6 

13,67 

846 

160,5 

13,71 

889 

157,5 


30 

8 

10,35 

25060 

332,7 

8,19 

4083 

111,7 

8,19 

3348 

96,0 

8,24 

4593 

127,6 


30 

10 

6,75 

18130 

179,8 

5,49 

4232 

77,7 

5,49 

3426 

65,9 

5,49 

3.64 

76,5 


Average 

9,40 

5450 

66,8 

6,10 

1263 

41,9 

6,10 

1203 

40,9 

6,10 

1231 

43,2 


15 

3 

17,18 

556 

2,0 

15,28 

467 

2,9 

15,28 

526 

3,3 

15,28 

510 

2,9 


15 

5 

8,36 

705 

1,5 

8,02 

562 

1,9 

8,02 

432 

1,7 

8,02 

441 

1,8 


15 

8 

6,82 

736 

1,0 

6,71 

669 

1,3 

6,71 

686 

1,2 

6,71 

733 

1,3 


20 

3 

17,00 

1370 

10,9 

14,55 

737 

15,6 

14,55 

1027 

15,9 

14,55 

1008 

16,6 


20 

8 

7,54 

2709 

7,7 

7,31 

1876 

7,7 

7,31 

1684 

7,1 

7,31 

2041 

7,9 


20 

10 

4,28 

2574 

5,7 

4,21 

2168 

6,2 

4,21 

2367 

6,4 

4,21 

2720 

7,1 

1 

25 

3 

16,95 

1167 

28,5 

15,41 

1016 

47,5 

15,41 

1055 

50,6 

15,41 

800 

50,9 


25 

8 

7,74 

5038 

30,4 

7,03 

1903 

22,9 

7,03 

2066 

23,7 

7,03 

2664 

26,7 


25 

10 

9,76 

13830 

66,2 

9,69 

13360 

71,2 

9,69 

10670 

57,1 

9,69 

13710 

70,3 


28 

3 

17,10 

7531 

168,8 

14,25 

2024 

112,8 

14,25 

1943 

114,0 

14,25 

2380 

119,8 


28 

8 

10,63 

20850 

188,2 

10,34 

11240 

133,7 

10,34 

12110 

144,9 

10,34 

10960 

128,0 


28 

10 

9,11 

93700 

705,6 

8,89 

66920 

519,5 

8,89 

61180 

522,6 

8,89 

69100 

527,7 


30 

3 

21,47 

4873 

190,1 

18,10 

4082 

295,5 

18,10 

2724 

255,0 

18,10 

2702 

239,1 


30 

8 

14,36 

258800 

3026,0 

13,96 

135200 

1701,0 

13,96 

124200 

1576,0 

13,96 

123600 

1526,0 


30 

10 

10,50 

103700 

898,1 

10,26 

55810 

497,4 

10,26 

61050 

566,9 

10,33 

65660 

606,6 


Average 

11,90 

34543 

355,4 

10,90 

19869 

229,1 

10,90 

18915 

223,1 

10,90 

19935 

222,2 


15 

3 

16,19 

789 

1,5 

14,68 

647 

2,1 

14,56 

546 

2,3 

14,56 

437 

2,2 


15 

5 

8,39 

489 

0,7 

7,71 

525 

1,1 

7,71 

487 

1,0 

7,71 

513 

1,1 


15 

8 

7,85 

366 

0,3 

7,79 

342 

0,5 

7,79 

417 

0,6 

7,79 

384 

0,5 


20 

3 

17,60 

929 

6,8 

15,17 

978 

11,8 

15,02 

648 

11,7 

15,02 

861 

12,8 


20 

8 

9,01 

4476 

8,4 

8,68 

4139 

8,5 

8,59 

3603 

8,5 

8,59 

3434 

8,5 


20 

10 

5,95 

1444 

2,6 

5,83 

1302 

3,6 

5,80 

1251 

3,8 

5,80 

1464 

4,2 

S 

25 

3 

16,83 

1486 

22,5 

15,84 

1671 

39,2 

15,74 

1077 

38,9 

15,74 

1148 

39,4 

1 

25 

8 

8,40 

15390 

59,5 

8,24 

12770 

57,5 

8,19 

11340 

53,7 

8,19 

11400 

54,1 


25 

10 

12,14 

1251 

5,6 

12,07 

1382 

9,5 

12,07 

1249 

8,9 

12,07 

1155 

9,0 


28 

3 

17,48 

4458 

80,8 

15,77 

2371 

88,3 

15,60 

1727 

83,5 

15,60 

1721 

82,9 


28 

8 

9,97 

16910 

119,1 

9,63 

11070 

93,3 

9,39 

9072 

81,6 

9,39 

7802 

74,6 


28 

10 

8,92 

10490 

49,7 

8,76 

9234 

55,4 

8,76 

10090 

57,1 

8,76 

10920 

63,8 


30 

3 

22,05 

3925 

147,3 

19,80 

3477 

200,0 

19,69 

2634 

189,6 

19,69 

2233 

177,0 


30 

8 

14,99 

52040 

514,7 

14,54 

33800 

348,4 

14,52 

29400 

326,1 

14,52 

29980 

339,6 


30 

10 

9,17 

36480 

263,3 

8,63 

26750 

205,3 

8,63 

26950 

229,9 

8,63 

27860 

228,9 


Average 

12,30 

10062 

85,5 

11,50 

7364 

75,0 

11,50 

6699 

73,1 

11,50 

6754 

73,3 
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Table 6: Computational resu’ 


is for FormulationTu^, 


with Valid Inequalities (II) bis 



n 


Pl2 

+V.I.j28).{35)-l36) 

e 1 2 + V . I . . (HD - HD 





RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 

RGAP 

Nodes 

time 


15 

3 

16,90 

701 

1,1 

16,32 

574 

1,6 

16,32 

576 

1,6 

16,32 

513 

1,6 


15 

5 

8,47 

759 

1,0 

8,45 

769 

1,3 

8,45 

850 

1,4 

8,45 

850 

1,5 


15 

8 

8,83 

556 

0,6 

8,81 

577 

0,8 

8,78 

582 

0,9 

8,78 

437 

0,9 


20 

3 

16,70 

876 

4,5 

15,94 

579 

7,1 

15,94 

649 

7,0 

15,94 

705 

6,9 


20 

8 

11,58 

4408 

7,9 

11,57 

4061 

7,9 

11,56 

3902 

8,3 

11,56 

3725 

7,8 


20 

10 

8,82 

2044 

4,7 

8,82 

2044 

5,4 

8,82 

2739 

6,2 

8,82 

2389 

5,6 

1 

25 

3 

14,56 

983 

13,6 

14,27 

689 

19,0 

14,27 

669 

19,8 

14,27 

750 

19,9 

1 

25 

8 

11,33 

30540 

107,2 

11,29 

29000 

103,6 

11,29 

25290 

96,5 

11,29 

28560 

104,6 


25 

10 

16,16 

5271 

19,2 

16,15 

4339 

18,4 

16,15 

4335 

18,6 

16,15 

5226 

21,4 


28 

3 

17,55 

5190 

86,6 

16,87 

2660 

74,6 

16,87 

2788 

79,9 

16,87 

2997 

83,0 


28 

8 

15,77 

25740 

160,5 

15,75 

26910 

175,6 

15,75 

25320 

179,3 

15,75 

26240 

173,8 


28 

10 

14,53 

53050 

269,0 

14,52 

57390 

296,1 

14,52 

58690 

314,7 

14,52 

50760 

265,1 


30 

3 

19,73 

5994 

151,6 

18,36 

4419 

183,1 

18,36 

3236 

148,0 

18,36 

4846 

190,2 


30 

8 

14,63 

108300 

965,3 

14,59 

119500 

1068,0 

14,59 

110300 

1094,0 

14,59 

133400 

1228,0 


30 

10 

11,77 

94490 

675,7 

11,77 

94490 

678,4 

11,77 

91310 

700,3 

11,77 

107000 

806,9 


Average 

13,80 

22594 

164,6 

13,60 

23200 

176,1 

13,60 

22082 

178,4 

13,60 

24560 

194,5 


15 

3 

17,58 

783 

1,9 

14,72 

678 

3,3 

14,49 

502 

3,2 

14,49 

594 

3,4 


15 

5 

7,85 

6345 

1,2 

6,80 

437 

1,8 

6,58 

387 

1,8 

6,58 

358 

1,9 


15 

8 

4,34 

162 

0,4 

4,04 

166 

0,7 

4,04 

150 

0,8 

4,04 

153 

0,8 


20 

3 

18,22 

1777 

10,7 

15,52 

959 

13,3 

15,25 

847 

14,1 

15,25 

904 

13,8 


20 

8 

7,09 

2821 

6,7 

6,77 

1272 

6,1 

6,62 

1260 

6,7 

6,62 

1001 

6,4 


20 

10 

4,03 

3045 

5,4 

3,84 

2437 

6,0 

3,74 

1879 

5,6 

3,74 

1950 

5,5 

I 

25 

3 

17,47 

1639 

28,0 

16,17 

1404 

46,7 

16,11 

1035 

45,3 

16,11 

1377 

50,7 

u 

u 

25 

8 

7,29 

3264 

23,5 

6,46 

1895 

23,2 

6,30 

1777 

23,8 

6,30 

1571 

23,1 

« 

25 

10 

7,28 

11330 

48,8 

7,09 

6548 

38,5 

6,99 

5529 

36,0 

6,99 

5793 

37,2 


28 

3 

17,16 

4146 

108,7 

14,52 

1557 

83,6 

14,21 

1076 

89,4 

14,21 

904 

82,9 


28 

8 

9,39 

14120 

139,8 

8,345 

5935 

89,7 

8,14 

5073 

82,3 

8,14 

5264 

84,9 


28 

10 

7,54 

34770 

247,2 

6,66 

26570 

206,3 

6,51 

11410 

113,9 

6,51 

9299 

99,6 


30 

3 

21,72 

4437 

194,9 

17,34 

3598 

269,3 

17,02 

2889 

237,2 

17,02 

3989 

279,3 


30 

8 

12,14 

114200 

1343,0 

10,75 

17250 

291,0 

10,66 

97400 

1434,0 

10,77 

16090 

303,5 


30 

10 

7,54 

30780 

305,4 

6,60 

13560 

170,7 

6,55 

8198 

120,1 

6,59 

8883 

133,7 


Average 

11,10 

15194 

164,4 

9,70 

5618 

83,3 

9,50 

9294 

147,6 

9,60 

3875 

75,1 


15 

3 

15,86 

2248 

4,4 

13,58 

2925 

6,9 

13,55 

2611 

6,7 

13,55 

2388 

6,7 


15 

5 

7,98 

877 

2,7 

7,35 

785 

3,5 

7,35 

898 

3,7 

7,35 

762 

3,4 


15 

8 

6,23 

726 

1,1 

6,14 

678 

1,4 

6,11 

804 

1,6 

6,11 

698 

1,5 


20 

3 

17,41 

4007 

21,0 

14,06 

4384 

31,1 

14,06 

4642 

31,9 

14,06 

4712 

32,7 


20 

8 

7,40 

5217 

15,8 

7,12 

4646 

15,9 

7,11 

5109 

17,8 

7,11 

4380 

15,7 


20 

10 

3,98 

3505 

7,8 

3,91 

2882 

8,5 

3,89 

3016 

8,8 

3,89 

3319 

9,5 

CO 

25 

3 

16,28 

4434 

63,9 

14,65 

5324 

98,9 

14,65 

4923 

91,3 

14,65 

6696 

102,0 

CJ 

o 

nJ 

oa 

25 

8 

8,07 

7836 

64,9 

7,39 

4695 

51,8 

7,39 

5066 

54,5 

7,39 

5400 

57,9 

25 

10 

9,99 

35520 

158,8 

9,96 

34470 

151,7 

9,96 

49260 

222,5 

9,96 

37030 

169,3 


28 

3 

17,11 

5934 

210,7 

14,05 

4436 

229,0 

14,05 

5004 

290,9 

14,05 

3642 

206,0 


28 

8 

9,14 

43450 

580,5 

8,61 

54300 

1787,0 

8,61 

35590 

559,7 

8,61 

32250 

503,8 


28 

10 

8,61 

140900 

1323,0 

8,20 

94410 

1013,0 

8,20 

105800 

1150,0 

8,20 

125300 

1328,0 


30 

3 

22,20 

9661 

517,0 

18,45 

7648 

623,9 

18,45 

7766 

637,6 

18,45 

6555 

539,6 


30 

8 

12,74 

187400 

4854, 

12,26 

177000 

4807, 

12,17 

176000 

4400, 

12,19 

184600 

4494, 


30 

10 

9,30 

253000 

3586, 

9,06 

191400 

2931,0^^' 

9,04 

208800 

3333,0<^^ 

9,04 

194700 

2991,0^'-' 


Average 

11,50 

46981 

760,8 

10,30 

39332 

784,0 

10,30 

41019 

720,7 

10,30 

40829 

697,4 

TOTAL 

11,70 

22470 

266,2 

10,40 

16108 

231,6 

10,30 

16536 

230,7 

10,30 

16197 

217,6 






















